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摘 要 : 替代 稳 态 下 的 荒漠 生态 系统 生态 弹性 反映 了 系统 承受 环境 干扰 后 的 恢复 能 力 ,这 对 认识 荒漠 生态 系统 过 
程 具有 重要 的 理论 意义 。 当 前 研究 主要 集中 在 湖泊 和 森林 生态 系统 ,而 有 关 苇 漠 生 态 系统 弹性 的 研究 普遍 未 考虑 
替代 稳 态 ,也 鲜 有 考虑 弹性 的 时 间 变 化 。 本 研究 以 量 康 北部 的 古 尔 班 通 古 特 沙漠 南 缘 至 沙漠 腹地 样 带 为 例 ,运用 
2001 年 1 月 至 2020 年 12 月 的 MODIS 全 球 植被 指数 遥感 数据 ,采用 BFAST(Breaks for Additive Season and Trend ) il 
状态 空间 建 模 对 数据 进行 处 理 和 提取 ， en linda 出 不 同时 段 奉 代 稳 态 条 件 下 的 生态 弹性 ， 
同时 根据 计算 结果 分 析 了 该 地 区 生态 弹性 时 空 演变 特征 及 其 影响 机 制 。 研 究 结果 表明 :(1) 研究 区 弹性 :和 弹性 整 
体 呈 现 先 下 降 后 上 升 的 趋势 ， gn 间 差 异 显著 ,(2) 生态 弹性 对 降水 变化 存在 消 后 响应 。(3) 降水 
季节 变化 的 差异 会 降低 降水 量 与 生态 弹性 之 间 的 相关 性 。 综 上 所 述 ,生态 弹 性 的 空间 分 布 总 体 受 降水 格局 控制 ， 
但 立地 条 件 导致 的 植被 空间 异 质 性 增加 了 生态 弹性 空间 分 布 的 复杂 性 ,而 生态 弹性 与 降水 变化 的 关系 取决 于 植被 
群落 构成 .植物 对 降水 变化 响应 .降水 量变 化 趋势 和 季节 分 布 。 本 研究 在 认识 荒漠 生态 系统 功能 稳定 性 维持 机 制 
以 及 荒漠 生态 保护 和 修复 方面 具有 重要 的 理论 与 实践 意义 
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生态 弹性 是 衡量 生态 系统 稳定 性 的 重要 指标 ， ”表意 义 。 研 究 表明 ,该 地 区 植被 对 气候 变化 反应 十 
其 定义 为 在 转换 到 一 个 由 不 同 结构 和 功能 的 替代 ”分 敏感 “, 过 去 数 十 年 呈现 恢复 -退化 -恢复 交替 变 
稳 态 之 前 ,系统 所 能 承受 干扰 的 大 小 '“。 以 系统 变 ”化 趋势 ”, 意 味 着 该 区 域 在 不 同 替代 稳 态 间 反 复 转 
量变 化 趋势 为 基准 ,系统 变量 围绕 变化 趋势 上 下 波 ”化 。 在 变化 环境 中 ,生态 弹性 是 生态 系统 维持 功能 
动 , 其 实际 数值 偏离 变化 趋势 的 时 间 越 长 ,生态 弹 稳定 的 重要 机 制 。 以 植被 功能 的 时 间 序 列 为 切入 
性 越 弱 。 而 系统 变量 存在 多 种 偏离 变化 趋势 的 情 点 ,探讨 旱 康 北 部 荒 江 生态 系统 不 同 替 代 稳 态 下 生 
Ob ,并 会 表现 出 稳定 性 ,这 就 是 蔡 代 稳 态 。 系 统 变 。 态 弹 性 的 时 空 属 性 特征 以 及 生态 弹性 内 在 的 作用 
量 受 正面 环境 扰动 时 高 于 变化 趋势 变化 的 替代 稳 机制, 在 荒漠 生 态 系统 过 程 认 识 ` 区 域 生 态 保 护 和 
ee 反之 为 稳定 状态 11”“。 不 同 站 ”修复 方面 兼 具 重 要 的 理论 和 实践 意义 。 


代 稳 态 下 生态 系统 表现 出 的 弹性 不 同 ,其 系统 变量 地 球 大 数据 相关 技术 的 研究 和 构建 ,为 利用 遥 
e 明显 ,这 种 现象 普遍 存在 于 湖 RAN a Pe EAS ARSC EEE T REE. 
泊 、 和 森林 荒漠 等 各 类 生态 系统 ”。 在 区 分 固有 的 季节 变化 和 环境 干扰 方面 最 重要 的 方 


古 尔 班 通 十 特 沙漠 是 发 育 典型 的 温带 车 漠 系 法 就 是 BFAST(Breaks for Additive Season and Trend) 方 
统 ,沙漠 南 缘 植 被 受 一 定 的 绿洲 排水 补给 ,沙漠 腹 IES. OTTER 导 到 
地 植被 主要 依赖 少量 降水 生存 。 量 康 北部 的 荒漠 ”广泛 应 用 ,并 可 用 于 量化 生态 系统 弹性 """。 通 过 时 
植被 ,其 空间 格局 在 整个 古 尔 班 通 古 特 沙 演 具有 代 间 序 列 量化 生态 系统 弹性 的 代表 性 方法 分 为 两 关 
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一 类 是 自 相 关系 数 法 , 另 一 类 是 自 回 归 系 数 法 。 前 者 
通过 计算 系统 变量 的 时 间 自 相关 系数 量化 弹性 后; 
后 者 则 是 计算 系统 变量 自 回归 系数 量化 弹性 …” 。 
两 种 方法 均 建 立 在 生态 系统 只 存在 单一 稳 态 的 前 
提 下 ,所 以 无 法 量化 奉 代 稳 态 下 的 生态 弹性 。 最 
近 , 生 态 系统 从 不 同 稳定 状态 退出 所 用 的 时 间 解 决 
了 替代 稳 态 条 件 下 生态 弹性 的 量化 问题 ,并 应 用 于 
湖泊 生态 系统 藻类 生物 量 的 生态 弹性 量化 中 ”， 
但 是 尚未 应 用 于 空间 尺度 上 的 生态 弹性 研究 。 
Meng 等 利用 Landsat 时 间 序 列 量化 出 不 同时 期 湘 
江 流 域 生态 弹 性 ,并 对 弹性 的 时 空 变化 做 了 初步 探 
索 , 但 是 缺少 对 替代 稳 态 下 弹性 的 差异 性 和 其 变化 
驱动 因子 的 深入 探讨 。 

为 了 深入 认识 旱 康 北部 荒 江 生态 系统 弹性 的 
时 空 演变 ,本 文选 择 了 一 条 从 沙漠 南 缘 至 沙漠 腹地 
的 样 带 ,运用 2001 年 1 月 至 2020 年 12 月 MODIS 的 
全 球 植被 指数 遥感 数据 产品 ,将 该 时 间 序 列 产品 分 
成 四 段 不 同时 期 ,并 使 用 BFAST 和 状态 空间 建 模 对 
数据 进行 处 理 和 提取 ,通过 退出 时 间 法 计算 了 替代 
稳 态 条 件 下 的 生态 弹性 ,并 进一步 分 析 降 水 变化 与 
生态 弹性 时 空 演 变 的 内 在 机 制 。 研 究 能 够 为 量化 
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1 研究 区 概况 与 方法 


11 研究 区 概况 及 样 带 选择 

研究 样 带 南 起 旱 康 北部 古 尔 班 通 古 特 沙 江南 
缘 , 北 至 沙漠 腹地 ,地 理 位 置 位 于 87°42'~88°12'E、 
44°21'~45°21'N( 图 1) ,南北 长 111.3 km, 东西 宽 
40.1 km, 总 面积 约 4.47x10’km*。 样 带 内 发 育 南北 
走向 的 固定 半 固 定 沙丘 ,沙丘 相对 高 度 10~30 m, FÉ 
带 内 沙丘 与 丘 间 低地 相间 而 列 ,地 表 起 伏 。 样 带 内 
降水 由 南 至 北 梯度 递减 ,区 域 年 平均 降水 量 在 100~ 
250 mm 2” 。 春 季 积 雪 融 水 是 土壤 水 的 重要 补给 2 ， 
也 是 惹 漠 植 物 的 重要 水 源 ”。 区 域内 发 育 以 梭 梭 
(Haloxylon ammodendron ) ARI (H. persicum ) 等 为 
建 群 种 的 典型 温带 荒漠 植被 ,常见 灌木 包括 淡 枝 沙 
HÆ (Calligonum leucocladum ) 等 是 该 区 沙丘 呈 固 定 
半 固 定 状 态 的 主要 因素 ”。 淮 木 以 外 ,区 域内 草本 
层 植物 稀 玻 分 布 , 稼 见 一 年 生 草 本 植物 如 猪 毛 菜 
(Solsola spp.) FARZ ( Ceratocarpus arenarius)55。 样 


图 1 研究 区 位 置 及 概况 


Fig. 1 Location and overview of the study area 
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带 南部 .中 部 .北部 生物 多 样 性 指数 呈 递 减 趋势 ” ， 
旦 荒漠 植物 物种 组 成 不 同 ,因此 沿 南北 方向 对 样 带 
进行 三 等 分 (图 la) ,并 进一步 分 析 降 水 与 生态 弹性 
时 空 演 变 。 
1.2 研究 方法 
1.2.1 生态 弹性 量化 ”生态 弹性 的 量化 流程 如 图 2 
所 示 。 文 中 数据 来 源 于 GEE(Google Earth Engine ) , 
调用 的 遥感 数据 为 MODIS MOD13Q1 级 NDVI 产 品 
数据 ,时 间 分 辨 率 为 16 d, 空 间 分 辨 率 为 230 m, 地 
FHL AL ER ZN WGS84, 由 于 产品 数据 始 于 2000 年 2 
月 ,为 保证 时 间 序 列 数据 完整 性 , 故 研 究 时 段 为 
2001 年 1 月 一 2020 年 12 月 。 本 研究 利用 最 大 值 合 
成 法 对 NDVI 数 据 进行 校正 ,通过 提取 同一 月 份 不 
同时 间 节 点 上 最 大 的 NDVI 数 值 ,作为 该 时 期 的 ND- 
VI 值 ,从 而 进一步 消除 环境 干扰 ,最 后 得 到 共 20 
a NDVI 的 月 时 间 序 列 数据 。 

生态 弹性 作为 长 期 观测 得 到 的 稳定 性 指标 , 需 
要 较 长 的 时 间 序 列 进行 量化 。 该 区 域 植被 指数 最 
高 值 出 现在 2016 年 ,图 2 显示 2015 年 前 后 到 2020 
年 NDVI 趋势 变化 较 前 期 具有 较 大 差异 , 故 将 
2016 一 2020 年 数据 划分 为 一 个 时 间 段 ,并 利用 该 时 
段 数 据 量化 同时 期 弹性 。 为 比较 不 同时 期 弹性 的 
变化 ,研究 中 量化 所 需 的 时 间 序 列 数据 需 等 长 , 结 
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2000 2005 2010 2015 2020 
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合 地 球 表层 系统 分 析 方 法 ””, 故 以 5 a 为 一 个 周期 进 
行 划 分 ,将 数据 以 5 a 为 一 个 周期 ,分 为 wt、ts 4 
个 时 间 段 ,分 别 代 表 2001 一 2005 年 .2006 一 2010 年 、 
2011—2015 年 .2016 一 2020 年 。 其 中 , 样 带 内 每 一 
个 顶 格 具体 的 时 间 序 列 为 也 ,其 中 也 由 下 列 公 式 中 
多 项 构成 : 

Y,=S,+T,+E, (1) 
式 中 :5 为 季节 项 ,表示 NDVI 时 间 序 列 随 季 节 谈 化 
的 函数 项 ;7 为 趋势 项 ,表示 NDVI 因 季节 变动 外 的 
压力 扰动 所 产生 变化 的 函数 项 ;为 残余 项 ,表示 
NDVI 因 脉冲 扰动 所 产生 变化 的 函数 项 ,本 人 研究 使 
用 BFAST 方 法 提取 以 上 各 项 。 

将 获取 的 NDVI 数 据 转 化 成 时 间 序 列 对 象 后 ， 
数据 处 理 的 主要 步 又 如 图 2 所 示 。 首 先 ,提取 不 同 
时 间 段 NDVI 时 间 序 列 ,利用 BFAST 检 测 变化 趋势 ， 
其 中 的 季节 项 使 用 默认 的 调和 模式 来 拟 合 。 之 后 ， 
利用 趋势 项 和 残余 项 进行 状态 空间 建 模 。 该 建 模 
有 效 地 反映 自然 界 的 真实 状态 ” ,具体 建 模 方法 调 
用 R 语 言 相 应 的 MARSS 包 2 。 完 成 状态 空间 建 模 
后 ,再 次 使 用 BFAST 分 离 趋势 项 和 残余 项 ,并 通过 
残余 项 划分 出 稳定 状态 I1(E,>0) 和 稳定 状态 II(E< 
0)。 最 终 ,利用 残余 项 计算 退出 时 间 ,从 而 得 到 生 
态 弹 性 。 


WOR 
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图 2 数据 处 理 流 程 及 残余 项 与 替代 稳 态 间 的 关系 


Fig.2 Methods for data processing and relationship between remainder component and alternative stable states 


202305.00244v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


5} PAE SE BFAS FLEE TIRER SHE FE 25 Je 811 
退出 时 间 计 算 的 是 环境 扰动 使 生态 系统 状态 ”应 时 间 段 内 弹性 及 其 变化 趋势 做 相关 分 析 。 降 水 


退出 当前 稳定 状态 所 需 的 平均 时 间 ”“”。 由 于 状态 
空间 建 模 还 原 了 数据 变化 的 连续 性 ,这 使 研究 可 以 
准确 地 计算 出 系统 在 不 同 奉 代 稳 态 之 间 转 换 所 用 
的 “退出 时 间 (EXITTIME)”。 这 里 借助 统计 方法 将 
不 同时 段 的 EXITTIMBE; 作 为 样本 数据 ,然后 计算 各 
阶段 平均 值 ,方法 如 下 : 

2, month, 


EXITTIME, = >, (2) 
i=l 


SUH i=l, …, n,n 表示 在 研究 时 间 内 系统 变量 于 同 
一 稳 态 内 做 连续 变动 的 频数 ,month; 为 第 i 次 连续 运 
动 内 所 包含 的 具体 月 份 数量 ,7 为 具体 的 时 间 段 。 
当 month; 中 <0 时 ,说 明 在 该 时 间 段 内 生态 系统 处 
于 NDVI 低 于 预期 的 稳定 状态 工 ,计算 结果 为 
EXITTIME,” ,弹性 -为 UEXITTIME;  ” ; 当 
>0 时 ,说 明 在 该 时 间 段 内 生态 系统 处 于 NDVI 高 
于 预期 的 稳定 状态 1, 计 算 结果 为 EXITTIME%“”， 
弹性 为 VEXITTIME"”…“”“。 计 算得 出 的 EXITTI- 
ME: 越 长 , NDVI 恢复 到 预期 数值 的 时 间 越 长 ,生态 
弹性 越 小 。 

1.2.2 生态 弹性 变化 趋势 及 降水 因素 相关 分 析 方 
法 ”为 深入 认识 弹性 与 降水 波动 变化 的 关系 ,本 研 
究 选 取 了 不 同时 期 多 年 平均 年 降水 量 与 4 sto ts tal 
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数据 使 用 的 是 CHIRPS(Climate Hazards Group Infra- 
Red Precipitation with Station data) 日 降水 量 数据 ,其 
空间 分 辨 率 为 0.05° ,并 通过 GEE 获取。 经 分 区 统计 
提取 生态 弹性 不 同 分 区 平均 值 ,然后 利用 规则 网 格 
(图 1a) 提 取 多 年 平均 年 降水 量 数据 以 及 生态 弹性 
数据 并 计算 Spearman 相关 系数 中 。 和 后 态 弹 性 变化 
趋势 主要 通过 后 期 数据 减 去 前 期 数据 得 到 ,如 到 
时 间 段 弹性 变化 趋势 等 于 时 间 段 弹性 减 去 .时间 
段 弹性 ,降水 量 的 变化 趋势 同 理 。 


2 结果 与 分 析 


21 生态 弹性 时 空 分 布 

弹性 数值 及 其 变化 趋势 的 时 空格 局 如 图 3。 弹 
性 : 取 值 范围 在 0.075~0.500。 时段 弹性 ' 整 体 相 对 
较 高 ,4 时 段 弹 性 ' 整 体 相对 较 低 , 可 能 在 该 区 域 不 同 
时 间 段 降水 量变 化 趋势 影响 了 区 域 弹性 。 弹 性 * 高 
值 在 样 带 南部 分 布 较为 集中 ,这 主要 受降 水 量 由 南 
至 北 梯度 递减 影响 。 弹 性 - 取 值 范围 在 0.085~ 
0.417。& 时 段 弹 性 -数值 高 于 其 他 时 间 段 ,ty 到 时 间 
段 高 值 呈 斑 块 状 混杂 分 布 , 说 明 弹 性 的 变化 受 局 部 
生境 条 件 影响 较 大 。 弹 性 变化 趋势 显示 ,到 时 间 
段 内 样 带 南 部 弹性 ' 旺 大 幅 下 降 趋势 ,ws 到 4 时 间 有 段 


G) tAE 0) 2 弹性 (k) -4 弹性 + 
i ae ge od Bi À 


D ibh E (m) iH (n) tt HE 


弹性 变化 趋势 M 0 15km 
-0.5 0.5 


图 3 不 同时 间 段 弹性 数值 ( 左 ) 及 变化 趋势 ( 右 ) 的 时 空格 局 


Fig. 3 Spatial distribution of resilience value and trends in different periods 
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存在 大 幅度 上 升 趋势 。 不 同时 期 弹性 与 弹性 -变化 
趋势 均 存 在 较 大 差异 ,可 能 与 区 域 植 被 对 不 同类 型 
环境 扰动 存在 适应 性 差异 有 关 。 

综合 来 看 ,弹性 对 不 同类 型 环境 扰动 的 响应 在 
空间 分 布 和 时 间 变 化 上 均 存 在 差异 ,这 种 差异 可 能 
受降 水 量变 化 趋势 .局 部 生境 条 件 和 植被 对 环境 扰 
动 的 适应 性 共同 影响 。 
2.2 生态 弹性 与 降水 量 关 系 的 空间 差异 

样 带 整体 多 年 平均 年 降水 量 在 .到 4 时间 段 内 
出 现 大 幅度 上 升 。 将 样 带 从 南 到 北 等 分 成 3 个 部 分 
(图 La) ,图 4 为 样 带 不 同 区 域 各 时 段 的 平均 年 降水 
量变 化 , 样 带 南 部 和 中 部 呈现 先 升 高 后 下 降 的 趋 
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势 ,而 北部 则 保持 上 升 趋势 。 

整体 上 ,弹性 与 降水 量 之 间 具 备 显著 的 相关 性 
( 表 1)。 弹 性 除 与 本 时 段 内 降水 相关 外 ,也 与 前 期 
的 降水 显著 相关 。 如 样 带 整 体 在 时间 段 的 生态 弹 
性 ' 与 4 时 间 段 降水 量具 备 较 高 的 相关 性 ,相关 系数 
为 0.592。 但 这 种 相关 性 存在 局 部 差异 ,人 研究 基于 降 
水 数据 分 辩 率 大 小 构建 矢量 网 格 (图 la) ,经 规则 采 
样 各 部 提取 60 个 样 点 的 不 同时 期 降水 及 生态 弹性 
数据 平均 值 ,最 终 计 算得 出 降水 与 弹性 相关 系 的 时 
空 变 化 。 样 带 内 南部 区 域 与 降水 量 保持 显著 正 相 
关 , 但 样 带 中 部 及 北部 到 时 间 段 弹性 与 降水 量 存 
在 显著 负 相关 。 此 外 ,不 同 蔡 代 稳 态 条 件 下 的 弹性 
与 降水 量 相关 性 差异 较 大 。 如 样 带 南 部 2 时 间 段 的 


(b) 南部 降水 量 


多 年 平均 年 降水 量 /mm 


多 年 平均 年 降水 量 /mm 


图 4 样 带 年 平均 降水 量变 化 


Fig. 4 Variations in averaged annual precipitation 


表 1 弹性 与 降水 量 相关 性 


Tab.1 Correlation between resilience and precipitation 


整体 南部 中 部 北部 
Spearman 相关 系数 - 
同期 降水 ”前 期 降水 同期 降水 ”前 期 降水 同期 降水 ”前 期 降水 同期 降水 ”前 期 降水 
弹性 * -0.279 *" 一 0.107 一 -0.049 一 -0.464"" 一 
弹性 - -0.071 0.596” = -0.511 E -0.409 z 
弹性 0.604 ** 0.592" 0.741" 0.701 0.061 -0.086 -0.288" -0.084 
É 弹性 0.075 0.084" 0.336" 0.318 -0.010 -0.143 -0.511” -0.392™ 
弹性 0.018” -0.182 0.689 0.682 -0.373” -0.238” 0.016 0.145 
弹性 - 0.033” 0.190™ 0.197 0.183 -0.215" -0.252 -0.162 0.054 
弹性 0.281™ 0.387™ 0.192 0.244" 0.232 0.017 0.178 0.254" 
E 弹性 - 0.105 0.059 0.336” 0.322” 0.327" 0.320 0.015 -0.119 


注 :*** 表 示 P<0.001, ** 9278 P<0.01 ,* 9275 P<0.05, FIA. 
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弹性 ”与 不 同时 期 降水 量 之 间 的 相关 性 为 0.701~ 
0.741 ,而 弹性 与 降水 量 之 间 的 相关 性 则 低 于 0.34。 
样 带 整体 弹性 变化 趋势 与 降水 量变 化 趋势 的 
相关 性 较 低 ( 表 2) ,而 样 带 整体 弹性 变化 趋势 相 比 
局 部 与 降水 量变 化 趋势 之 间 存 在 更 显著 的 相关 
性 。 缘 合 来 看 ,弹性 "和 弹性 -与 降水 量 相关 性 随时 
间 和 空间 变化 复杂 多 样 ,可 能 与 区 域 植 被 对 降水 变 
化 趋势 存在 适应 性 差异 有 关 ,并 且 这 可 能 还 与 降水 
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季节 分 布 的 时 空 差异 有 关 。 
2.3 生态 弹性 与 降水 季节 分 布 的 关系 

图 5a 为 稳定 状态 1 和 稳定 状态 工 在 各 月 份 出 现 
的 频数 ,二 者 的 频数 差 反 映 了 整个 样 带 在 不 同月 份 
下 稳定 状态 变化 的 规律 。 整 个 生态 系统 有 更 大 概 
率 在 6 月 和 9 一 11 月 处 于 稳定 状态 1, 在 3 一 5 月 及 7 
月 处 于 稳定 状态 II( 图 5b)。 

当 降 水 量变 化 促进 植被 生长 时 (图 5c), 生 态 系 


表 2 弹性 与 降水 量变 化 趋势 的 相关 性 


Tab.2 Correlation between the trends of resilience and precipitation 


; 相关 系数 整体 南部 中 部 北部 
Spearman AR p PF = SEs = RA = = 
i 同期 降水 ”前 期 降水 同期 降水 ”前 期 降水 同期 降水 ”前 期 降水 同期 降水 前 期 降水 
弹性 * 0.070 一 -0.024 一 0.302" 一 -0.384" - 
hrt: 
弹性 - -0.122 一 0.137 一 -0.081 一 -0.155 一 
弹性 0.418" -0.027 0.171 -0.158 0.010 0.074 —0.288° 0.499" 
tots 
弹性 - -0.071 0.285" 0.156 0.133 0.213" 0.274" -0.326" 0.208 
弹性 * -0.260™ -0.387 -0.277° -0.231° 0.422" 0.316" -0.211' -0.027 
ta—ta 
弹性 - -0.111 0.350" -0.166 -0.089 -0.146 0.291% 0.139 0.246 
12.0 (a) 不 同 稳定 状态 出 现 频数 
as 0 稳定 状态 工 a BERET 
110 3 
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图 5 稳定 状态 频数 与 降水 量变 化 的 关系 


Fig .5 Stable state frequency in relation to precipitation variability 
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统 在 5 一 6 月 的 稳定 状态 1 所 停留 的 时 间 增 加 ,7 月 
在 稳定 状态 I 停 留 的 时 间 减 少 ,甚至 存在 转换 为 替 
代 稳 态 的 可 能 。 因 此 ,生态 系统 在 稳定 状态 I 的 平 
均 退 出 时 间 增 加 ,弹性 数值 降低 。 反 之 ,弹性 -数值 
降低 。 

当 春 夏 两 季 降 水 量变 化 对 稳定 状态 的 影响 有 
差异 时 (图 5e, 图 5f) ,稳定 状态 III 的 平均 退出 时 间 
将 更 加 复杂 ,与 不 同时 期 降水 量 的 相关 性 降低 。 
此 ,这 可 能 是 蔡 代 稳 态 下 的 弹性 与 部 分 时 段 的 降水 
量 之 间 相关 性 不 高 的 主要 原因 。 


3 讨论 


本 文 以 息 康 北部 荒漠 植被 样 带 为 例 , 通 过 结合 
BFAST 状态 空间 模型 和 退出 时 间 法 ,初步 量化 并 探 
讨 了 不 同 殖 代 稳 态 下 息 康 北部 荒漠 生态 系统 南 毕 
至 腹地 生态 弹性 的 时 空 变化 。 人 研究 表明 ,该 区 域 生 
态 弹 性 总 体 呈 先 下 降 后 上 升 的 趋势 ,不 同 替 代 稳 态 
条 件 下 弹性 数值 从 古 尔 班 通 古 特 沙漠 边缘 至 沙漠 
腹地 的 空间 分 布 差异 显著 , 且 弹 性 数值 对 降水 波动 
存在 滞后 响应 ,降水 波动 的 空间 差异 .植物 群落 的 
异 质 性 分 布 及 其 对 降水 变化 的 差异 响应 可 能 是 造 
成 这 种 差异 的 主要 成 因 。 

通过 NDVI 计 算得 出 的 殖 代 稳 态 下 的 生态 弹 
性 ,表明 了 生态 系统 功能 在 承受 正面 或 负面 的 环境 
扰动 后 恢复 的 时 间 长 度 。 生 态 系统 恢复 时 间 越 长 ， 
表明 其 弹性 越 小 ,承受 扰动 和 恢复 的 能 力 越 低 。 对 
芒 漠 生态 系统 而 言 ,降水 是 生态 系统 生产 力 的 最 主 
要 限制 因素 。 根 据 Bai 等 3 对 整个 中 亚 地 区 生态 系 
统 弹 性 研究 发 现 , 插 康 北 部 的 古 尔 班 通 古 特 荒漠 生 
态 系统 南 缘 生 态 弹性 大 于 沙漠 腹地 ,与 样 带 北部 弹 
性 :在 和 时 间 段 低 于 南部 的 研究 结果 一 致 。 

本 研究 结果 显示 沙漠 南 缘 至 腹地 降水 变化 趋 
势 差 异 明显 , 样 带 南 部 及 中 部 降水 量 呈 现 先 升 高 后 
下 降 的 趋势 (图 4b ,图 4c) ,北部 降水 量 则 一 直 保 持 
上 升 趋势 (图 4d) ,与 新 疆 大 部 分 地 区 一 致 拉 。 而 降 
水 量变 化 趋势 的 空间 分 布 差异 是 不 同 替代 稳 态 条 
件 下 的 弹性 在 tt 时 间 段 内 空间 分 布 差异 的 主要 
原因 ,结果 显示 弹性 变化 趋势 与 降水 量变 化 趋势 间 
显著 相关 。 如 样 带 整 体 弹性 ' 在 到 时 间 有 段 的 变化 
趋势 与 同时 间 段 降水 量 的 相关 性 为 0.418 ,这 说 明 弹 
性 会 随 着 同期 降水 量 的 上 升 而 上 升 。 

弹性 指标 反映 了 系统 在 历经 变化 时 保持 其 原 


有 功能 .结构 .特性 和 抗 干 扰 及 重组 的 能 力 。 因 而 , 生 
态 系统 内 物种 的 丰富 度 将 有 助 于 增加 生态 弹性 ” 。 
已 有 研究 表明 ,荒漠 植被 群落 的 弹性 受 物 种 丰富 度 
的 影响 显著 后 ,而 降水 量 高 的 区 域 植 物 物 种 丰富 度 
明显 大 于 降水 量 低 的 区 域 扩 。 降 水 由 沙漠 南 缘 向 
腹地 的 复杂 变化 导致 的 植物 物种 丰 度 变化 "3, 这 是 
息 康 北部 古 尔 班 通 古 特 沙漠 南 缘 生 态 弹 性 高 于 腹 
地 的 自然 基础 。 此 外 ,人 研究 表明 立地 条 件 如 坡 向 和 
坡度 也 会 对 生态 弹性 产生 影响 ”…。 乍 康 北 部 荒漠 
生态 系统 ,沙丘 纵横 起 伏 ,植物 稀 玻 分 布 , 且 主 要 分 
布 于 沙丘 底部 和 丘 间 地 i。 这 是 因为 地 形 的 起 伏 
导致 丘 间 地 有 更 好 的 水 分 条 件 , 使 丘 间 低地 具备 更 
高 的 物种 丰富 度 ,进而 影响 弹性 的 空间 分 布 ,而 不 
同类 型 沙丘 (如 固定 沙丘 、 半 固定 沙丘 与 流动 沙丘 ) 
由 于 水 分 条 件 差异 也 会 使 得 荡 漠 植被 种 群 的 稳定 
状态 在 不 同 地 貌 类 型 条 件 下 存在 差异 ~”。 因 此 ,在 
降水 格局 和 局 地 立地 条 件 综合 作用 下 ,里 康 北部 荒 
漠 生态 系统 弹性 的 空间 分 布 既 有 南北 方向 上 的 趋 
势 性 差异 ,在 局 部 又 呈现 弹性 斑 块 混杂 分 布 的 特征 
(图 3)。 

以 往 研 究 表明 古 尔 班 通 古 特 沙漠 降水 量 与 ND- 
VI 之 间 存 在 滞后 响应 Y, 本 文 研究 进一步 证 明 这 种 
滞后 响应 也 会 反映 到 基于 NDVI 计 算 的 生态 弹性 变 
化 。 如 样 带 南部 存在 弹性 与 前 期 降水 量 之 间 的 相 
关 性 大 于 0.7。 这 是 因为 降水 波动 对 同期 植被 物种 
丰富 度 的 影响 过 程 较 为 缓慢 ,通过 改变 植被 群落 结 
构 从 而 反映 在 生态 弹性 上 的 时 间 较 长 ,因此 生态 弹 
性 对 降水 量变 化 存在 滞后 响应 。 此 外 ,并 非 所 有 的 
荒漠 植物 盖 度 与 降水 量 呈 显著 正 相 关 52 ,如 和 白 梭 
TBS ` 淡 枝 沙 拐 束 等 的 密度 和 盖 度 随 降水 的 增加 而 减 
少 ,而 梭 梭 与 降水 则 没有 相关 性 ,这 可 能 与 植被 在 
适应 温带 荒漠 气候 环境 过 程 中 形成 了 不 同 生 存 策 
略 有 关 “' 。 本 文 研 究 结果 显 示 该 区 域 存 在 弹性 与 
降水 量 呈 负 相 关 的 现象 ( 表 1)。 如 样 带 北 部 在 六 到 
5 时 间 段 弹性 与 同期 降水 量变 化 趋势 之 间 的 相关 系 
数 在 -0.409 ~ -0.511。 结 合 弹性 与 降水 季节 分 布 的 
关系 可 知 , 当 样 带 中 部 及 北部 植被 的 盖 度 与 降水 量 
呈 负 相关 时 ,弹性 随 降水 量 增加 而 下 降 。 由 于 降水 
的 季节 分 布 变化 多 样 , 且 对 弹性 的 影响 过 程 复杂 , 
研究 区 荒漠 生态 系统 弹性 与 降水 量 的 相关 性 也 呈 
RASH. RAMA ,研究 区 生态 弹性 与 降水 量 的 
相关 性 取决 于 植被 群落 构成 、 降 水 的 变化 趋势 及 降 
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水 的 季节 分 布 。 

为 了 解决 不 同 替代 稳 态 条 件 下 的 生态 弹性 量 
化 问题 ,本 研究 采用 了 BFAST 变 ee 
SE 状态 (稳定 状 
态 I 和 稳定 状态 卫 ) ,并 利用 退出 时 间 法 量化 得 出 生 
态 弹 性 ,通过 息 康 北部 古 尔 班 通 古 特 荒漠 生态 系统 
研究 的 应 用 ,证 实 了 还 漠 生 态 系 统 普遍 存在 和 蔡 代 稳 
态 ,以 及 不 同 替代 稳 态 条 件 下 的 生态 弹性 的 时 空 
差异 。 与 已 有 人 研究 工作 ”3 相 比 ,本 文 研究 通过 划 
分 不 同 替 代 稳 定 状 态 和 不 同时 段 研究 荒漠 生态 系 
统 弹 性 ,更 好 地 解释 了 其 S 
化 规律 ,为 认识 荒漠 生 性 长 期 动态 与 环境 干 
扰 的 关系 提供 了 新 的 思路 。 阜 康 北 部 的 古 尔 班 通 
古 特 沙漠 稳定 性 是 保障 区 域 生 态 环境 安全 、 促 进 区 
域 可 持续 发 展 的 重要 因素 中 ,其 生态 弹性 的 研究 结 
果 有 助 于 深入 认识 环境 干扰 下 切 漠 生态 系统 功能 
ee Ey Ol AS a Te a 
复 也 具有 重要 的 理论 指 


70 态 稳定 


4 结论 


本 研究 运用 BFAST 和 状态 空间 建 模 的 方法 , 通 

过 对 遥感 时 间 序 列 数据 的 分 期 处 理 提 取出 了 不 同 
时 段 殖 代 稳 态 条 件 下 的 弹性 ,并 量化 不 同时 段 阜 康 
北部 的 荒漠 生态 系统 弹性 的 空间 分 布 ,然后 分 区 探 
讨 了 不 同时 段 生 态 弹 性 与 降水 变化 之 间 的 关系 。 
弹性 分 布 总 体 受 降水 格局 控制 ,但 立地 条 件 导致 的 
植被 空间 异 质 性 增加 了 弹性 空间 分 布 的 复杂 性 ;在 
时 间 上 ,生态 弹性 在 研究 时 段 内 时 先 下 降 后 上 升 的 
趋势 ,对 降水 变化 的 响应 存在 滞后 性 ,生态 弹性 与 
降水 量 的 相关 性 取决 于 植被 群落 构成 .植物 对 降水 
变化 的 响应 以 及 降水 量变 化 趋势 和 季节 分 布 。 本 
文 为 量化 奉 代 稳 态 条 件 下 的 生态 弹性 提供 了 方法 
与 技术 支撑 ,在 深入 认识 荒漠 生态 系统 功能 稳定 性 
维持 机 制 以 及 荒漠 生态 系统 保护 和 修复 方面 具有 
理论 与 实践 意义 。 
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Spatial and temporal patterns of ecological resilience under alternative stable 
states in the desert of the north Fukang region 


Xu Yuzhe'’, Lin Tao, Li Jun” 
(1. State Key Laboratory of Desert and Oasis Ecology, Xinjiang Institute of Ecology and Geography, Chinese Academy of 
Science, Urumqi 830011, Xinjiang, China; 2. University of Chinese Academy of Sciences, Beijing 100049, China; 
3. Technology Innovation Center for Desert -Oasis Ecological Monitoring and Restoration, Ministry of Natural Resources, 
Urumqi 830022, Xinjiang, China; 4. Akesu National Station of Observation and Research for Oasis Agro-ecosystem, 
Xinjiang Institute of Ecology and Geography, Chinese Academy of Science, Urumqi 830011, Xinjiang, China) 


Abstract: The ecological resilience of desert ecosystems under alternative stable states reflects the ecosystem’ s 
ability to recover from environmental disturbances, which has important theoretical implications for 
understanding the desert ecosystem processes. Recent research has focused on some intensively studied 
ecosystems such as lakes and forests, whereas studies related to the ecological resilience of desert ecosystems 
have hardly considered alternative stable states as well as the temporal changes in resilience. In this study, based 
on the MODIS global vegetation index from January 2001 to December 2020, we used Breaks for Additive 
Season and Trend and state-space modeling to determine the temporal- spatial pattern of ecological resilience 
under alternative stable states quantified by the exit times for a transect extending from the southern edge to the 
hinterland of Gurbantunggut Desert in north Fukang. Moreover, mechanisms underlying the temporal- spatial 
pattern of ecological resilience were explored. The results of this study indicated the following: (1) Resilience* 
and Resilience” showed a decreasing trend within the study period but with significant differences between the 
southern edge and the hinterland of the desert. (2) Ecological resilience also showed a lagged response to 
precipitation changes. (3) Differences in the seasonal patterns of precipitation could weaken the correlation 
between precipitation and ecological resilience. In summary, the spatial distribution of ecological resilience is 
generally controlled by the precipitation patterns, but the spatial heterogeneity of vegetation due to site conditions 
increases the complexity of the spatial distribution of ecological resilience, with the relationship between 
ecological resilience and precipitation resulting from the combined effects of the composition of vegetation, 
adaptation of plants to changes in precipitation, and the seasonal pattern of precipitation. Our study is important 
for understanding the mechanisms that desert ecosystems utilize to maintain functional stability as well as in 
desert conservation and ecological restoration. 


Keywords: alternative stable states; ecological resilience; exit time; Gurbantunggut Desert 


